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Abstract 

We have performed a theoretical study of the effects of pinning potential and dis- 
sipation on vortex tunneling in superconductors. Analytical results are obtained in 
various limits relevant to experiment. In general we have found that pinning and dissi- 
pation tend to suppress the effect of the vortex velocity dependent part of the Magnus 
force on vortex tunneling. 
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To address the question of the resistance of a superconductor at zero temperature, a clear 
understanding of vortex tunneling is needed. Despite ample experimental evidence for vortex 
tunneling in superconductors H, the theoretical study is of dimensional analysis in nature|| 
[3| and controversial j|] . Furthermore, in vortex tunneling there is no comprehensive study of 
the role played by the Magnus force, a force with a part proportional to the magnitude of 
vortex velocity but perpendicular to the direction of the velocity This is a very unsatisfactory 
situation, in view of the fact that the Magnus force is a general property of a vortex || and 
tunneling is extremely sensitive to dynamics. In this paper we present an analytical study 
on vortex tunneling in the presence of the Magnus force, with a special attention paid to the 
effects of pinning and dissipation. A general picture of those effects on vortex tunneling can 
be obtained from our study. 

We start with the Hamiltonian for a vortex. To make the presentation simple, we consider 
vortex tunneling in a semi-infinite superconductor film which defines the x — y plane. The 
supercurrent is along the x-direction, and the edge of the film is lying on the x-axis. The 
length scale in the present study is always much larger than the size of a vortex core. A 
vortex can then be regarded as a point particle, and as discussed in Ref. its effective 
Hamiltonian can be generally written as 
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with the vector potential A determined by V x A = hp s dz/2. The meaning of each term 
in eq.(l) is as follows. The vector potential reflects the existence of the the vortex velocity 
dependent part (VVDP) of the Magnus force F m = q v hp s d{y s — r) x z/2. The Magnus force 
depends on the relative velocity between the superfluid velocity v s and the vortex velocity r, 
and the superfluid velocity dependent part (SFVDP) of the Magnus force will contribute to 
the vortex potential V. In accordance with the calculation of tunneling, the vector potential 
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will be taken as 



A= -p s d (y, 0,0) 



(2) 



and it can be shown that results are independent of the choice of gauge. Here q v = +1(— 1) 
stands for the vorticity parallel (anti-parallel) to the unit vector z in ^-direction, h is the 
Planck constant, p s is the superfluid electron number density with the factor 1/2 counting 
for the Cooper pairing, and d is the thickness of the film. We take the vortex mass M to 
be finite as demonstrated in Ref.||, and will show that it is relevant to vortex tunneling. 
The vortex potential V contains both the contributions from the SFVDP Magnus force and 
pinning centers. In the following we shall take V to be of the form 



The pinning potential in the x-direction is approximated by the harmonic potential, and 
k x should be determined experimentally. The potential V\{y) consists of the contributions 
from the SFVDP Magnus force and the pinning potential in the ^/-direction, which has a 
metastable point at y = in the present paper. In the limit of no pinning and extreme type 



only the contribution from the SFVDP Magnus: the potential due to external supercurrent 
Vo and the image potential from the edge of the superconductor film. Here m c is the mass of 
a Cooper pair. The dissipative environment of the vortex, the last term in eq.(l), consists 
of a set of harmonic oscillators as formulated in Ref. 0. The jth harmonic oscillator has a 
mass rrij and a frequency Uj. The effect of the dissipative environment is specified by the 
spectral function 



V(r) = V x {y) + h x x 2 . 



(3) 
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In the present paper, we shall assume the spectral function to have the following form 

J(uS) = f]uj s exp < i , (6) 

I u c ) 

with uj c the cutoff frequency whenever needed. In accordance with Ref.0, s > 1 is the 
superohmic case, s = 1 the ohmic case, and < s < 1 the subohmic case. In the ohmic 
damping case, r\ is the vortex friction such as discussed in Ref. ||. 

A remark on the relationship between the present study and others. We note that the 
vortex motion is identical to the motion of an electron in the presence of a magnetic field 0. 
Results obtained in the present paper can then be directly applied to that case with q v 
replaced by electron charge e and hp s /2 by magnetic field 5[|TI|. For the case of a vortex in 
superfluid, the superfluid electron pair number density p s /2 should be replaced by superfluid 



helium atom number density p s [ 1 1 



The tunneling is described by the Euclidean action 0, ||, [TTJ] 



S= dr 
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where j3 = l/k B T is the inverse temperature. The tunneling rate is equal to exp{— S c /h}, 
where the semiclassical action S c is determined by the bounce solution of the equation 5S = 0, 
in which the periodic boundary condition, (r(h/3), {qj(fof3)}) = (r(0), {^-(0)}), is required. 

We are interested in the vortex tunneling out of the metastable state y = 0. After the 
tunneling in y-direction other degrees of freedoms can take arbitrary value. Therefore the 
summation over the final states, integrations over environmental degrees of freedoms, {qj}, 
and over the x coordinate, will be taken. Those are gaussian integrals. After the integrations, 
the effective Euclidean action is 



rh/3 M -11 r hf3 rh/3 

/ dr -Mf + V x {y) +-/ dr dr' [k(\r - r'|) + g(r - r')][y(r) - y{r')f 
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with the normal damping kernel k as 
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k(r) = -i duJ{u) . L ; * " , (9) 



and the damping kernel due to the VVDP Magnus force coupled x-direction motion, which 
we shall call the anomalous damping kernel, as 



g(r) = ~— i-P d\ f Mul + in (10) 



Here 
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and 



u l = 77 " ( 13 ) 
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In the large r limit, from eqs.(6) and (9) the normal damping kernel takes the form 

Kr) = - V 4n > ( 14 ) 

which demonstrates the one-to-one correspondence between the low frequence part of the 
spectral function and the long time behavior of the damping kernel. Now we have obtained 
an effective one-dimensional problem. This is the starting point for the following detailed 
study of the tunneling in various experimentally accessible limits. 

no pinning, no dissipation Without disorder and dissipation, there is no quantum decay 



of supercurrent in an infinite filmJT2[. In reality such as discussed in the present paper the 
film always has an edge, tunneling is possible. To see the role played by the VVDP Magnus 
force, we first study the case with no dissipation and no pinning. Because of the VVDP 
Magnus force, motions in the x- and y- directions become coupled to each other. Therefore 
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although the normal damping kernel k vanishes now, the n = mode in the anomalous 

kernel g contributes to the effective action. Then the effective action is 

2 „ 
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We note that in the language of spectral function, this corresponds to the subohmic bath 
case with s = 0(c./. eqs.(8) and (14)), and eq.(15) is explicitly gauge invariant under the 
change y — > y + constant, because of the periodic boundary condition of x imposed in the 
tunneling calculation. 

For low enough density p s , that is, for a weak VVDP Magnus force, there is a tunneling 
solution. Particularly, for very low density and low temperatures the semiclassical action 
may be evaluated perturbatively: 

2 
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where y c {r) is the bounce solution at zero temperature without the VVDP Magnus force. 
Eq.(16) shows a remarkable p 2 s dependence and linear temperature dependence [13 1. 

For high enough density, that is, for a strong VVDP Magnus force, the tunneling rate 
vanishes at zero temperature, because the VVDP Magnus force renormalizes the original 
potential such that the state near y = is stable. This can be noted from eq.(15) where as 
T — > the cross term in the last term disappears. A straightforward calculation leads to the 
criterion for the localization as 
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Thus we have obtained that the s = dissipative environment is marginal for localization 



in tunneling decay, compared to the s — 1 case for the tunneling splitting! 14]. This is the 



dynamical localization caused by the VVDP Magnus force |12|. It is worthwhile to point out 



that according to eq.(17) although a large VVDP Magnus force inhibits vortex tunneling, a 
large vortex mass instead favors vortex tunneling. 

no pinning, finite dissipation Now we set the pinning potential k x /M = oj 2 = in 
eq.(10). Both the normal damping kernel and the anomalous damping kernel are finite. 
Since the quantum tunneling corresponds to the limit of large imaginary time, H(3 — > oo, 
we look for the large time limit behavior of the anomalous damping kernel g. In the this 
limit, we may replace the summation l/hf3J2n by the integration 1/271"/ dv. Then for the 
environment < s < 2, we find the anomalous damping kernel g in the large r limit as 

9{T) = a 2^M^ Psd ) ^7^' (18) 
with a a numerical constant of order of unity. Here 

O n poo — \ 
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and J(u>) = rju s has been used. The effective dissipative environment corresponding to the 
anomalous damping kernel is s e ff = 2 — s(c.f. eqs.(6) and (14)), which leaves the ohmic 
damping unchanged, transforms the subohmic damping into superohmic damping and vice 
versa. 

For the case s > 2, using eq.(4) for the spectral function J, we find that the effective 
dissipative environment corresponding to the anomalous damping kernel is s e ff = 0, which 
smoothly connects the result for < s < 2. 

An important example is the ohmic damping case, where an exact expression for a in 
eq.(18) can be obtained. Carrying out a detailed but straightforward calculation, we find 
the effect action at zero temperature as 

/oo r 1 "| | roo roo | 
dr -my 2 + V^y) + — Veff / dr / dr' — [y(r) - y{r)f , (20) 
-oo lZ \ Z7T J-oo J-oo t — t r 
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with the effective damping strength T] e ff as 

= + (21) 

The effect of a dissipative environment can now be summarized as follows. It has been 



demonstrated in Ref. fllj] that the subohmic dissipation has strong effects on tunneling, while 
the superohmic dissipation has weak effects. From the above analysis we have that if the 
normal damping kernel k(r) of eq.(9) is subohmic, the anomalous damping kernel g(r) of 
eq.(10) is superohmic, and vice versa. Then according to Ref.[|14]] the effect of the VVDP 
Magnus force, represented by the anomalous damping kernel, is weak/strong on vortex tun- 
neling for the normal subohmic/superohmic damping. In particular, for normal superohmic 
damping with s > 2, the vortex tunnels as if there were no effect of dissipation. For the 
normal ohmic damping with s — 1, we need to compare the relative strength of the VVDP 
Magnus force and the dissipation according to eq.(21). Classically, the VVDP Magnus force 
tends to keep a vortex moving along an equal potential contour, but the friction instead 
along the gradient of the potential. In general, we can conclude that the dissipation tends 
to suppress the effect of the VVDP Magnus force on vortex tunneling. 

finite pinning, no dissipation The normal damping kernel vanishes in this situation. 
The anomalous damping kernel can be expressed by hyperbolic functions. We find that the 
effective action is then 

rh/3 M "I l ( h \ 2 rhf3 rhf3 

S,„ = J dr [-Me + VM\ + — -p,d\ j dr / dr'x 
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It is the superohmic case with s e ff = oo, because the effective spectral function has the form 
5(uj — u) x ), which has no low frequency mode. Therefore the tunneling rate is nonzero for any 
magnitude of the VVDP Magnus force. This result shows that pinning has a very strong 



influence on the vortex tunneling in the presence of the VVDP Magnus force, because the 
introducing of the pinning potential bends the straight line trajectory of a vortex and makes 
the transition to other trajectories possible. We note that by letting uj x = we recover 
eq.(15), where the tunneling rate vanishes for a sufficiently strong VVDP Magnus force. 

We can evaluate the semiclassical action perturbatively, if the density is low, or, pinning 
potential is strong as done in eq.(16). In the large density limit semiclassical action may be 
evaluated by a variational method similar to Ref . 0] in the case of ohmic damping. However, 
we have a more powerful method to perform the calculation, because the strong VVDP 
Magnus force freezes the kinetic energy of the vortex. In this case x and y coordinates now 
form a pair of canonically conjugate variables, and the Euclidean action is||, 



rh/3 

S eff = j Q dr 



i^psdxy + Vi{y) + ^k x x 2 



(23) 



Following the calculation outlined in Refs. |], [K| , we find the semiclassical action in the form 



S c = hp s d J^dyJ^^-, (24) 



with y t the turning point determined by the equation V\{y) = PRlgl- The result shows that 
as k x increases, the semiclassical action decreases. Therefore we conclude that the pinning 
potential in the x-direction helps the tunneling in y-direction[[J. 

general case: finite pinning and dissipation In general we need to go back to eqs.(8-10) 
to study the tunneling under the influence of pinning and dissipation. However, based on 
the insight gained by the above analysis here we can draw a general conclusion without a 
detailed study: Since both the dissipation and pinning tend to suppress the effect of the 
VVDP Magnus force, the total effect of them will do the same. Indeed the anomalous 
damping kernel is always superohmic. In particular, in the presence of strong pinning and 
ohmic damping, the superohmic-like anomalous damping kernel g may be ignored compared 
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to the ohmic normal damping kernel k. Then from eq.(8) we have the effective action as 

r h/3 m in r h P r n P 1 

S eff = j o ^[-M^ + V^J+^jf drj o dr'j-—^[y(r)-y(r')]\ (25) 

which looks as if there were no effect of the VVDP Magnus force. This may explain the 
pronounced experimental observation of the absence of the effect of the VVDP Magnus force 
in vortex tunneling experiments on dirty superconductors 0. 

To summarize, we have performed a complete analytical study of the influence of pin- 
ning and dissipation on vortex tunneling. The VVDP Magnus force tends to decrease the 
tunneling rate. In the absence of pinning and dissipation, there is no tunneling for a strong 
VVDP Magnus force. Detailed results have been obtained when one or all of them is absent. 
Both pinning and dissipation tend to suppress the influence of the VVDP Magnus force on 
vortex tunneling. Present results may explain the absence of effects of the VVDP Magnus 
force in vortex tunneling experiments on dirty superconductors. 
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